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Abstract: Mitsui and Ohshima (2008) criticized the power-stroke model for muscle 
contraction and proposed a new model. In the new model, about 41% of the myosin heads 
are bound to actin filaments, and each bound head forms a complex MA3 with three actin 
molecules Al, A2 and A3 forming the crossbridge. The complex translates along the actin 
filament cooperating with each other. The new model well explained the experimental data 
on the steady filament sliding. As an extension of the study, the isometric tension transient 
and isotonic velocity transient are investigated. Statistical ensemble of crossbridges is 
introduced, and variation of the binding probability of myosin head to Al is considered. 
When the binding probability to Al is zero, the Hill-type force-velocity relation is resulted 
in. When the binding probability to Al becomes finite, the deviation from the Hill-type 
force-velocity relation takes place, as observed by Edman (1988). The characteristics of the 
isometric tension transient observed by Ford, Huxley and Simmons (1977) and of the 
isotonic velocity transient observed by Civan and Podolsky (1966) are theoretically 
reproduced. Ratios of the extensibility are estimated as 0.22 for the crossbridge, 0.26 for 
the myosin filament and 0.52 for the actin filament, in consistency with the values 
determined by X-ray diffraction by Wakabayashi et al. (1994). 
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1. Introduction 

In 1999, Mitsui [1] criticized the power-stroke model on the muscle contraction mechanism and 
proposed a new model. In 2008, Mitsui and Ohshima [2] refined the new model and discussed the 
steady filament sliding in detail demonstrating that the calculation results were in good agreement with 
experimental observations. They also outlined the discussion on the isometric tension transient and the 
isotonic velocity transient given in [1] by citing some calculation results. Thereafter, however, a few 
readers of [1] commented that it was very difficult to understand the discussion on the transient 
phenomena in [1], since there was no detailed explanation on the molecular processes upon which the 
theoretical treatment was based. In the present paper, we have largely revised that part of [1], trying to 
make the discussion more readable. 

Now the article [2] is regarded as Part I of Remarks series, in which the basic ideas of our model 
are introduced and steady muscle behaviors are discussed. In the present paper (Remarks II), 
non-steady muscle behaviors are discussed. We are preparing an article as Remarks III, in which 
discussion will be done on more recent experimental studies as cited in the articles [3,4]. 

The basic ideas of our model introduced in [2] are summarized as follows. A simple thermodynamic 
relation is derived, which indicates that there is an inconsistency in the power stroke model or 
swinging lever model. Our model is proposed to avoid this difficulty. It is assumed that a myosin head 
forms a complex with three actin molecules when it attaches to an actin filament. Here it should be 
noted that Andreeva et al. [5] found the evidence that the crossbridge can interact with more than 
1 actin monomer. The complex corresponds to the crossbridge. According to the X-ray diffraction 
studies [6-9], the intensity ratio of the [1, 0] and [1,1] equatorial reflections increases only minimally 
as the shortening velocity increases, indicating that the total number of myosin heads in the vicinity of 
the actin filament decreases only slightly. Taking this fact into account, it is assumed that about 41% 
of the myosin heads forms the crossbridges at any sliding velocity. Then mutual cooperativity takes 
place among the crossbridges in filament sliding, so that energy dissipation becomes reasonable 
magnitude (the order of kT) for one step of the crossbridge movement. Calculation based upon the 
model well reproduce the force-velocity relation given by Hill [10] and the energy liberation rate vs. 
force relation given by Hill [11]. 

Since the present study is based upon ideas that are quite different from the power stroke model and 
others, the basic ideas of our study are explained in some detail in Section 2. In Section 3, discussion is 
done of how variation in the crossbridge binding affects the muscle tension. In Section 4, discussion is 
done on the difference between the molecular processes in Phases 1 of the isometric tension transient 
and of the isotonic velocity transient, and extensibility ratios for the crossbridge, the myosin filament 
and the actin filament are estimated. Time course of the isometric tension transient is studied in 
Section 5. Time course of the isotonic velocity transient is studied in Section 6. The deviation from the 
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Hill-type force-velocity relation is derived in Section 7. Obtained results are summarized and 
discussed in Section 8. 

2. Basic Ideas for Discussion of the Transient Phenomena 

2.1. Deviation from the Hill-Type Force-Velocity Relation 

In the present study, the deviation from the hyperbolic force-velocity relation shown in Figure 1 has 
very important implication. The red line is our calculation results reported in [2] which agrees with the 
empirical hyperbolic force-velocity relation proposed by Hill [10]. The curve will be called Hill-type 
below. Edman [12], however, reported that carefully measured velocity deviated from the Hill-type 
relation as given by the black circles in Figure I, which he called a double-hyperbolic force-velocity 
relation. In Figure 1 , the tension T at which the experimental data start to deviate from the Hill-type 
curve is indicated as and the tension T at which T becomes 0 is denoted as T 0o \, s . The following 
values are determined by the data presented in Figure 6 A of [12]. 

T d JT 0 = 0.66 (1) 

WJo = 0.88 (2) 

Figure 1. Force-velocity relation and definition of Tdev and Tbobs. The red curve is the 
calculation result reported in [2], which coincides with the hyperbolic force-velocity 
relation proposed by Hill [10]. Circles are the experimental results cited from Figure 6A of 
the paper by Edman [12]. 
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2.2. MA3 Complex and U*n Transition 

In our model, about 41% of the myosin heads are bound to actin filaments (cf. Equation 3-1-1 
in [2]), and each bound head forms a complex MA3 with three actin molecules. The complex MA3 
translates along the actin filament changing the partner actin molecules. Figure 2 illustrates the step 
motion of MA3. The translation of the crossbridge along the actin filament is made possible by the 
structural change of MA 3 induced by the force fj which exerts on the junction J between the 
crossbridge and the actin filament. 
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Figure 2. MA3 complex (shadowed) and its step motion along the actin filament (after 
Figure 6 of [2]). A myosin head is expressed as neck domain plus catalytic domain and the 
actin molecules in MA3 as 1, 2 and 3. (a) Just before the myosin head moves to the right; 
(b) Just after the head (and the complex MA 3 ) moves to the new position; (c) The head is 
ready for next movement to the right. (Note that the actin molecule bound to a myosin 
molecule is called 2.) 
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Figure 3 shows the distribution of potential of force exerted on the myosin head in MA3. There are 
three potential wells on the actin filament corresponding to the three binding sites Al, A2 and A3. The 
wells are also as called Al, A2 and A3. Figure 3(a) shows the potential distribution considered in [2], 
where a myosin head exists solely in well A2. It is assumed that the potential barrier U* depends upon 
the force fj by U* = U*o - afi (Equation 3-5-2 in [2])). The kinetics of the myosin head will be 
discussed based upon Eyring's theory of the rate process (cf. [13]) as was done by Huxley and 
Simmons [14] and in [2]. Then the probability that the myosin head moves from well 2 to 3 across the 
potential barrier U* becomes proportional to exp(— U*lkT). Calculation was done in [2] based upon the 
model in Figure 3(a). An example of calculation results is cited by the red curve in Figure 1. Since the 
calculation well explains the experimental data for T < T^ v . this model seems close to reality when 
T < Tdev The deviation from the Hill-type relation when T > Td ev suggests that the model in Figure 3(a) 
should be modified, and the model in Figure 3(b) becomes the object to be considered in the present 
study. It is assumed that myosin heads start to move over the potential barrier U*n when T becomes 
larger than T^, and heads are distributed in Al and A2 when T > Td ev . The potential peak between Al 
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and A2 will be called U*n and the potential barrier for backward movement from A2 to Al will be 
denoted as U*\ 2 b and the barrier for forward movement from Al to A2 as U*nf- In the followings, 
transitions of these types will be called U*n transition. 

Figure 3. Potential of the force exerted on the myosin head in the MA 3 complex. The 
shape varies depending upon the tension T. (a) The case of T < Td ev , where the force-velocity 
relation is Hill-type. The black circle indicates that the myosin head exists solely in well A2; 
(b) The case of Tdev < T < T 0o b s , where the existence probability of the myosin head is 
distributed in wells Al and A2, causing deviation of the force-velocity relation from Hill-type. 
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The displacement of the myosin head from A2 to Al was briefly mentioned in Section 3.6 of [2] in 
term of "pull-up transition". In the article [1], the isometric tension transient was discussed in relation 
with the U*n transition by the statistical mechanics. Although the present study is based upon the 
same idea, we have found that such statistical-mechanical approach as in [1] makes the discussion very 
complex. We shall discuss the problem in a different manner as described below. 

2.3. Definition of Crossbridge Shortening Y 

In the following discussion, a parameter y is frequently used to represent shortening (which is 
positive for negative length change) of the crossbridge as in [2]. Unfortunately, the symbol y was used 
to represent elongation of the crossbridge in [1]. In the present paper, as an extension of the discussion 
of [2], y represents shortening. Figure 4 is to make the definition of y clear. K and J, respectively, are 
ends of the crossbridge on the myosin and actin filaments. The symbol x is defined as the projection of 
the vector from K to J on the myosin filament. Since the positions J are set on the right-handed red 
dotted line indicated as z in Figure 4, x is determined by the position of K. As shown in Figure 4(a), x 
is denoted as x eq when the tilting angle of the myosin head neck domain is at the equilibrium angle 
Equation- Then y is defined by 

y = x - x eq (3) 

The red dotted line indicated as x eq in Figure 4 indicates the origin of y. Thus y = 0 in Figure 4(a), 
y < 0 in (b) and y > 0 in (c). As discussed in [2], the crossbridge with shortening y exerts the force p(y) 
on the myosin filament. The stiffnesses for forward and backward forces are denoted as Kf and /q,, and, 
according to Equation 3-4-6 in [2], p(y) is given by 

p{y) = -Kfy for y < 0 (4a) 
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p(y) = -/<hy for y>0 



(4b) 



Values of aj and /c b are given in Appendix. 

Figure 4. Definition of y, the shortening of the crossbridge: y = x - x eq (after Figure 4 of [2]). 
The binding positions J are set on the dotted red line z. (a) The myosin head is at its 
equilibrium angle, 0 eq . x = x eq and y = 0; (b) The myosin head is pulling the myosin 
filament forward. The crossbridge is elongated and y < 0; (c) The myosin head is pushing 
the myosin filament backward. The crossbridge is shortened and y > 0. 
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2.4. Time Constant of the U* Transition 

To discuss the transition phenomena, it is important to have an idea about the mean time interval of 
occurrence of the U* transition. In the steady filament sliding, the mean time interval is equal to the 
mean time r s t ep in which a myosin head moves from one actin molecule to the neighboring one. This 
r step is given by 

r ste p = Llv (5) 

where L is the distance between the centers of actin molecules 2 and 3 in Figure 2 (the value of L is 
given in Appendix) and v is the sliding velocity. Values of r ste p are calculated by using v given by the 
red curve in Figure 1, and shown as a function of the relative tension in Figure 5, which indicates that 
r ste p is about 2 ms at T/T 0 = 0, becomes 10 ms around T/T 0 = 0.5 and increases to 1 s around T/T 0 = 1. 
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This result means that U* transition generally does not contribute to the early processes in the 
transition phenomena. 

Figure 5. Time interval Z" s t ep for a step motion of the myosin head as a function of T/Tq. 
r ste p = L/v, which is approximately equal to the time interval of the occurrence of the 
U* transition. 
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2.5. Statistical Ensemble of Crossbridges and Crossbridge Binding Probability p 

In discussion of the steady filament sliding [2], the crossbridge shortening y repeats the cycle from 
y c -L to y c through filament sliding and from y c to y c -L through the U* transition (cf, Figure 7 in [2]), 
where y c is the parameter which decreases as the tension T increases (as shown in Figure 8 given later.) 
In the following discussion, this temporal scheme is replaced by a spatial scheme. In our model, the 
ratio r = (number of myosin heads bound to actin filaments)/ (total number of myosin heads) is 
constant independent of tension (r = 0.41, cf. Appendix). Accordingly, we can consider an ensemble of 
bound crossbridges of definite number independent of the tension. 

Suppose that all complexes MA3 in right half sarcomeres are collected and superposed putting A2 
at the same position. This ensemble can be characterized by the ratio p(y)dy = (number of crossbridges 
having shrinkage between y and y + dy)/(total number of crossbridges). By definition, integration of 
piy) is 1 . Actually p(y) is a smooth function. The mean value of y in the p distribution is denoted as 
<y> and y c is defined by 

y c = <y> + (1/2) (6) 

Then p(y) is expected to be large between y c -L and y c . In the following discussion, p{y) is 
approximated to be constant (1/L) between y c -L and y c and zero outside of the region. This will be 
called rectangular p approximation. Equivalent approximation was used in the temporal scheme in [2] 
for studies of the steady filament sliding. 

The force-velocity relation deviates from Hill-type for T > 7d ev as shown in Figure 1 . The origin of 
the deviation is considered as follows. In Figure 3, the potential of force exerted on the myosin head in 
the MA3 complex is shown. The potential at A2 is lower than that at Al, and all the myosin heads are 
present in well A2 in the case of T < 7d ev as shown in Figure 3(a). On the other hand, as seen in 
Figure 4(b), the backward force exerted on the myosin head becomes stronger as y negatively increases, 
so that the myosin heads having negatively large y tend to bind to Al through the U*u transition. 
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Accordingly, the distribution of the myosin heads becomes as symbolically shown in Figure 3(b) for 
T> Tdev This effect is regarded as the origin of the deviation of force-velocity relation from the Hill-type. 

Symbols voobs and v c dev are defined as y c 's at T 0o \, s and Tdev, respectively, on the Hill-type 
force-velocity relation (cf Figure 8 given later.) As mentioned above, y c decreases with increasing T, 
and thus T < corresponds to y c > y c dev- 

Figure 6 shows examples of the rectangular p distribution for the steady filament sliding. 
Figure 6(a) is for fast sliding and (b) is for relatively slow sliding in the case of y c > y c d ev or T < Tdev 
Figure 6(c) is for slow sliding in the case of y c < y c d e v or T > Td ev . In Figures 6(al), (bl) and (cl), the 
upper horizontal line is an axis of the crossbridge shortening y and the black thick segment on it shows 
the regions of finite p(y) of length L. The lower horizontal lines correspond to the actin filament and 
Al, A2 and A3 represent actin molecules. The crossbridges having negative y produce positive stress 
p(y) according to Equation 4a. The red segments with their ends at A2 are examples of such 
crossbridges. The crossbridges having positive y produce negative stress p{y) according to Equation 4b. 
The blue segments with their ends at A2 or Al are examples of such crossbridges. 

Figure 6. Statistical ensemble of crossbridges in the steady filament sliding. 
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Figures 6(a2) and (b2) show the p distributions as functions of y, p(y) by the rectangular 
approximation, where p(y) is equal to (1/L) between y c -L and y c , and equal to zero outside of the 
region. Crossbridges in the red area produce positive stress, and crossbridges in the blue area produce 
negative stress. As the edge y c shifts to the left from (a2) to (b2), the red region increases and the 



Int. J. Mol. Sci. 2011, 12 



1705 



crossbridge ensemble produce more stress and the filament sliding becomes slow (Note tCf > /q, 
in Appendix). 

Figure 6(c) corresponds to the case of T > Td ev and thus y c < y c dev In (cl), the blue dotted segment 
symbolically indicates that a portion of crossbridges of small y bind to Al and produces negative 
tension, while the red dotted segment indicates that the rest of crossbridges having the same small y 
still bind to A2 and produce positive tension. In (c2), the violet triangle symbolically indicates that the 
crossbridges in this triangle are bound to Al. As shown in (c3), the violet triangle is replaced by the 
violet rectangle containing the same number of the crossbridges, to make calculation easy. In the 
following calculation, it is assumed that all the crossbridges in the violet rectangle are bound to Al and 
produce negative tension. 

Now let us calculate values of y c dev and y c oobs- By definition, y c d ev and y c oobs are y c 's for r 0o bs and Td ev 
in the Hill-type force-velocity relation. The mean tension per one crossbridge is denoted as p in [2], 
which can be obtained by integration of the tension p(y). Calculation is done by using p(y) = -K{ylL 
(Equation 4a) in the red region and piy) = ~n\,y/L (Equation 4b) in the blue region in Figure 6 (a2) or 
(b2). Result of the integration is 

p = {-(-/Cf/2)(y c - Lf + (-/r b /2)y c 2 }/L (7) 

This relation leads us to the Hill- type force-velocity relation in cooperation with the U* transition [2]. 
The isometric tension in the Hill-type relation is denoted as po and T 0 following [2]. Then, by using 
p/po = T/To (Equation3-l-2 in [2]), the relative tension T/T 0 is given by 

T/T 0 = {-(-Kf/2)(y c - Lf + (-/r b /2)y c 2 }/(Lp 0 ) (8) 

This equation gives relations between Tdev and y c dev and between r 0o bs and y c oobs- Combining 
Equations 1 , 2 and 8 gives 

ycdev = 1.60 nm, (9) 
ycoobs = 1 .03 nm. (10) 
In Figure 6, y c dev and y c o 0 b s are indicated by the vertical dotted red lines. 

3. Effect of U*n Transition on the Tension 

Equation 8 leads us to the Hill-type force-velocity relation. The stress T and y c related by 
Equation 8 is denoted as T 11 and y c H , which are in agreement with the observation when y c H > y c d ev . The 
observed T and y c are denoted as T* and y c * when y c H < y c d ev . 

Figure 7 illustrates the p distributions in various cases of y c H < y c dev Figure 7(al) shows the p 
distribution by Equation 8 having the edge y c H . This p distribution gives 7^ even though y c < y c d ev . It is 
an imaginary state in which the U*n transition is absent. Actually, however, the U*n transition takes 
place and the violet area appears causing decrease of the red area as shown in (a2), where the width of 
the violet area is denoted as x(y c ). Naturally T in (a2) is smaller than T in (al). Then y c in (a2) is 
changed into y c * in (a3) to make T the same as T in (al), where the width of the violet region is 
denoted as x. The change of the blue area from y c to yc* causes an increase of the tension while the 
change of the violet area from x(y c ) to x causes an decrease of the tension. If the effect of the former 
change is larger than that of the latter change, the tension in (a2) will increase and can be the same as 
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in (al). Below we shall assume that such tension adjustment actually takes place and discuss the 
relation between y c * and T*. 

Figure 7. Compensation of the effect of the U*n transition by shift of y c H to y*. 
Magnitudes of y c o 0 bs and y c dev are exaggerated for illustration, (a) y c o 0 bs < y c H < y c dev and 
y c * > 0; (b) y* = 0; (c) y c * < 0. 
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Since the tension is the same in Figures 7(al) and (a3), the blue area in (al) is the same as the sum 
of the blue area and violet area in (a3). Figures 7(bl)-(b3) show the case of y* = 0. Similarly to (al) 
and (a3), the blue area in (bl) is the same as the violet area in (b3). The width x of the violet area in 
(b3) is denoted as xo. Figure 7 (c) shows the case that y c * becomes negative. The brown area means 
that crossbridges in this area bind to Al and produce positive stress. The width of the violet area is 
indicated as x. Extrapolating the change from x in (a3) to xo in (b3), x in (c) is assumed to be larger 
than xo. Then sum of the red area and brown area in (c) is smaller than red area in (b3) and thus the 
tension produced in (c) is smaller than the tension in (b3). Then the stress produced in (b3) is the 
maximum that the muscle machine can produce, and should be equal to Tbobs- Accordingly, y c H in (bl) 
should be equal to v c oobs- Since the red area in (bl) is the same as the red area in (b3), y c B in (bl) is 
equal to xq in (b3). Thus we have 



Xo — VcOobs 



(11) 
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The tension T* deviates from T H at y c 



ycdev, and reaches its maximum at y c * = 0. To express 



such characteristics of T* vs. y* relation, the following set of equations are used: 

T*IT Q = -ay c * 2 + b 

a = (r 0o bs - 7dev)/(7oy c dev 2 ) = °- 086 O™ 1 ) -2 
b = T 0o bs/To= 0.$$ 



(12) 
(13) 
(14) 



In Figure 8, the blue curve illustrates the T*IT 0 vs. y c * relation given by this set of equations, and 
the red curve shows the Hill-type T^/To vs. y c u relation. The blue curve deviates from the red curve at 
v c = yc* = Jcdev and exhibits its maximum at y c * = 0 where T*/T 0 = 7o 0 bs/7b- 

Figure 8. Relative tensions as functions of y c . Red line: T^/To as a function of y c u . Blue 
line: T*IT$ as a function of y c *. (For definition of y c0 and y c (0), see Appendix). 
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Up to here, T/To is used for the relative tension, where To is the maximum tension in the extended 
Hill-type force-velocity relation. The experimentally observed maximum tension, however, is r 0o b s 
where r 0o bs = 0.88Tb (Equation 2). Hereafter discussion is concerned with experimental data and 
relative tension is defined as T/Toobs- In our model T/To = p/p 0 (Equation 3-1-2 in [2]), and the force p 
corresponding to r 0o b s is denoted as poobs- Then 



POobs _ 0.88 po 



(15) 



and 



?77oobs = p/poobs (16) 

If we put y c = ycOobs in Equation 7 and use xo = y c oobs (Equation 11), we get theoretical expression of 

Poobs as 



Poobs = {-(-Kf/2)(-L + xof + {-K h /2)x 0 2 } IL 



(17) 
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4. Phases 1 in the Isometric Tension Transient and Isotonic Velocity Transient, and Extensibility 
Ratios for the Crossbridge, Myosin Filament and Actin Filament 

Huxley [15] divided the transient responses to the sudden reduction of length or of load into four 
Phases. In our model, however, there is no exact correspondence between molecular processes of the 
four Phases in the tension transient and those in the velocity transient. To avoid confusion, we use the 
terms, Phase Tn in the tension transient and Phase Vn in the velocity transient, where n = 1, 2, 3, 4. 

Phase Tl is simultaneous decrease of tension in the isometric tension transient, while Phase VI is 
simultaneous shortening of muscle in the velocity transient. The length change per half sarcomere is 
denoted as AL\ ls in both transients. The experimentally measured isometric tension is Tbobs and the 
relative tension is defined by T/T 0obs . The relative load used in the paper by Civan and Podolsky [8] is 
expressed as the relative tension T/T 0o bs below. Figure 9 shows the experimental results on 
Phases Tl and VI, by circles for the isometric tension transient cited from the paper by Ford et al. [16] 
and by squares for the isotonic velocity transient cited from the paper by Civan and Podolsky [17]. 

Figure 9. T/T 0o \, s vs. ALh s relations in Phase Tl and Phase VI. Circles: experimental data 
cited from Figure 13 of Ford et al. [16]. Squares: experimental data obtained from Figure 3 
of Civan and Podolsky [17], in which the error bar means error of read-out from Figure 3 
of [17]. The green curve: r/r 0o bs calculated by Equations 18 and 20 with Ccbt = 2.2 
(Equation 22). The brown curve: T/T 0ohs calculated by Equations 19 and 21, 
with Ccbv = 4.6 (Equation 23). 
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At first sight, it was puzzling to see that distributions of the experimental data are quite different for 
the two cases, since the structural changes seem to be purely elastic both in Phases Tl and VI. Then it 
was reminded that the length change of sarcomere is a sum of those of the crossbridge, myosin 
filament and actin filament, which are proportional to each other (cf. the review by Irving [18]). This 
means that there are three elastic components, and it is a possibility that they have different response 
times in elastic changes from each other. 

In this connection, the experimental data reported by Julian and Sollins [19] seem important. They 
measured tension changes of single frog skeletal muscle fiber at increasing speed of step shortening. 
Their Figures 2 and 4 show experimental data on relative force vs. ALh s relation at different speeds of 
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shortening. The distribution of open circles in their Figure 2 (the length change period of about 1 ms) 
is similar to that of the data by Civan and Podolsky (1-2 ms) cited in our Figure 9. The distribution of 
filled triangles in Figure 4 of [9] (length change period of 0.4 ms) is close to the data by Ford et al. 
(the length change period of 0.2 ms) cited in our Figure 9. These facts seem to indicate that there are 
two kinds of elastic process in the tension response: The fast one occurs within about 0.4 ms and the 
slow one occurs between about 0.4 and 1 ms after the length change. The fast process seems 
responsible to the change in Phase Tl and combination of the fast and slow processes seems 
responsible to Phase VI. 

It seems plausible that elastic response of crossbridge and myosin filament almost simultaneously 
occurs since they belong to the same molecule and elastic response of actin filament occurs with some 
delay. Thus it is assumed that crossbridge and myosin filament are responsible to the change in Phase 
Tl and that all three components are responsible to the change in PhaseVl. Now changes of y c in 
Phases Tl and VI should be different from each other even when the length change ALh s is the same. 
Calculation is done on this assumption by using the rectangular p distribution shown in Figure 10. 

Figure 10. The p distributions before and after instantaneous changes in the isometric 
tension transient and in the isotonic velocity transient, (a) The isometric tetanus state; 
(b) Just after Phase Tl in the isometric tension transient; (c) Just after Phase VI in the 
isotonic velocity transient for the same length change as (b). 
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Figure 10(a) shows the p distribution at the isometric tension (the same as Figure 7 (b3)). 
Figure 10(b) shows the p distribution just after Phase Tl and (c) the one just after Phase VI for the 
same length change ALh s as (b). The violet areas in Figure 10 (b) and (c) are the same as (a) since U*n 
transition does not occur yet. Changes of the edge v c in Phases Tl and VI are denoted as Ay cT i and 
Avcvi- For the same ALh s , the tension T in Phase Tl is smaller than T in Phase VI in Figure 9. The edge 
Avcxi in Figure 10 (b) is set larger than Av c vi in (c) so as to make the red area in (b) smaller than that in 
(c), in accordance with the fact that the tension T in Phase Tl is smaller than T in Phase VI . 
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Since elastic changes of the crossbridge, myosin filament and actin filament are proportional to each 
other [10], the Ay cT i in (b) and Ay cV i in (c) are proportional to AL hs , and expressed by 

Ay cT i = -AWCcbt (18) 

and 

Aj c vi = -AWCcbv (19) 

where C C bt and Ccbv are constants. 

Let the tensions in Figure 10(b) and (c) be denoted as Tn and 7Vi, respectively. They are given by 
integration of -Kb{y - (-L)} in the violet area, -Kf_y in the red area and ~K\,y in the blue area. Thus we 
have 

7Vr 0o bs = [(-W2) {(x 0 + Avon) 2 - Avcxi 2 } - ("K f /2)(-L + *.oAy C Ti) 2 + (-K b /2)Ay cT i 2 ]/(Lpoobs) (20) 
and 

7Vi/r 0o bs = [(-K b /2){(x 0 + Ay cV i) 2 - Ay cV i 2 } - ("K f /2)(-L + x Q Ay cV i) 2 + (-K b /2)Ay c vi 2 ]/(Lp 0 obs) (21) 

Calculations were done for various trial values of Ccbt (Equation 18) and Ccbv (Equation 19). The best 
fit for the experimental data are obtained for the values 

Ccbt = 2.2 (22) 

Ccbv = 4.6 (23) 

Figure 9 shows the calculation results for r T i/7o 0 bs by the green curve and for 7Vi/7o 0 bs by the brown 
curve. They are in good agreement with the experimental data. 

Let us denote the ratios of extensibilities of the crossbridge, the myosin filament and the actin 
filament as tcb, r M , and r A at the elastic equilibrium respectively. Shortenings of the three per half 
sarcomere are denoted, respectively, as AycB, Ay M and Ay a, and the length change of half sarcomere as 
AL hs . Then, at the elastic equilibrium, they are given by 

AycB = -rcBAL hs (24a) 

Ay M = -r M AL hs (24b) 

Ay a " r A AL hs (24c) 

Naturally, 

rcB + r M + r A = 1 (25) 

On the above assumption, only the length changes of the crossbridge and the myosin filament 
contribute to Phase Tl. Then we have 

Ay cX i= ALh s / Ccbt = {tcb/^cb + r M )}AL hs (26) 

Thus, 

Ccbt = (rcB + r M )/rcB (27) 

As assumed above, the load change period is long enough and elastic changes of the crossbridge, 
myosin filament and actin filament contribute to Phase VI. Thus, Ay cV i is given by 
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Aj c vi = AL hs /CcBV = {rcEs/(/cB + r M + r A )}AL hs = r C BAL hs (28) 

Hence, 

Ccbv = 1/rcB (29) 
By using the values of Ccbt (Equation 22) and Ccbv (Equation 23), Equations 25, 27 and 29 give 

r CB = 0.22 (30a) 

r M = 0.26 (30b) 

r A = 0.52 (30c) 

The extensibility ratios were investigated by X-ray diffraction by Huxley et al. [20] and 
Wakabayashi et al. [21]. The values reported by Wakabayashi et al. [21] are tcb = 0.31, tm = 0.27 and 
r A = 0.42. Considering approximate nature of the theory and experimental errors, our values in 
Equation 30 are in reasonable agreement with the X-ray values. 

In this section, it is assumed that elastic change of the actin filament does not take place during the 
fast length change in Phase Tl. Then the elastic change of the actin filament should contribute to the 
next step, the tension recovery in Phase T2. This problem is discussed in next section. 

5. Isometric Tension Transient 

As noted in Section 4, Huxley [15] divided the transient responses to the sudden reduction of length 
into four phases. Time courses of these Phases are as follows. Phase 1 is instantaneous drop of tension. 
Phase 2 is rapid early tension recovery in next 1-2 ms. Phase 3 is extreme reduction or even reversal 
of rate of tension recovery during next 5-20 ms. Phase 4 is the gradual recovery of tension, with 
asymptotic approach to the isometric tension. 

As in Section 4, we use the terms Phase Tn (n = 1, 2, 3, 4) in the tension transient to avoid 
confusion. Phases Tn are more related with molecular process rather than time sequence. Phase Tl is 
simultaneous drop of tension caused by elastic shortening of the crossbridge and myosin filament. 
Phase T2 is due to elastic shortening of the actin filament. Phase T3 is related with the U*n transition, 
and divided into T3a and T3b. Phase T3a is the first part of Phase T3. Both Phases T2 and T3a 
contribute to the rapid early tension recovery (Phase 2 of Huxley). Phase T3b is the second part of T3 
where the extreme reduction or reversal of rate of tension recovery occurs (Phase 3 of Huxley). 
Phase T4 is the gradual recovery of tension, with asymptotic approach to isometric tension (Phase 4 of 
Huxley). (There were misprints in Section 4.4 of the article [2], and the last three sentences of the 
section should be neglected.) 

Figures 11 and 12 illustrate how the rectangular p distribution changes during these Phases. 
Figure 1 1 is for the case of Av c ti < y c de V and Figure 12 is for the case of Ay cT i > y c dev, where Ay cT i is y c 
just after Phase Tl. (cf. Figure 10(b)). 

As mentioned above, it is assumed that the elastic changes of the crossbridge and myosin filament 
occur in Phase Tl and then the elastic change of actin filament occurs in Phase T2. Figure 11(a) 
represents the states just after Phase Tl. Figure 11(b) shows the state just after Phase T2. There is a 
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time lag for the U*n transition to occur, and its effect is neglected in Phases Tl and T2, so that the 
width of the violet area is kept the same as the isometric value xo. 

Figure 11. Changes of the p distribution during the isometric tension transient in the case 
of Ay C Ti < Vcdev Magnitude of y c dev is exaggerated for illustration. 
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As noted referring to Figure 5, the mean time (r ste p) needed for the U* transition is relatively large 
(e.g., larger than 10 ms for T/Tq > 0.5). It is assumed that the potential barrier U*nt (£ forward) is 
lower than U* as shown in Figure 3 (b) and thus the U*uf transition starts before the U* transition. 
Phase T3a is regarded as the state that the U*nf transition is present but the U* transition does not 
occur yet, while the U*nt an d U* transitions coexist in Phase T3b. In Phase 4, only the U* transition 
exists. Figure 11(c) and (d), respectively, show the states at the end of Phases T3a and T3b. The violet 
area in (c) is smaller than in (b) because a portion of the crossbridges in the violet area change their 
binding partners from Al to A2 through U*nt transition. Thus the red area increases and the stress is 
stored and the tension increases in Phase T3a. The U* transition starts at (c) and the stress stored 
during Phase T3a is released by the shift of the edge y c from y cT 3a in (c) to y C T3b in (d). The shift of the 
edge from y cT3a (c) to y cT 3b (d) causes increase of the blue area and decrease of the violet area. These 
two effects tend to cancel each other and reduce change of the red area in the time course from (c) to 
(d), i.e., in Phase T3b. The reduced change of the red area may cause the "extreme reduction of rate of 
tension recovery" mentioned at the beginning of this section. In Phase 4, the U* transition and the 
filament sliding continue and the state approaches to the tetanus state (e). 

Figure 12 shows changes of the p distribution in the isometric tension transient for the case of 
AycTi > Vcdev Figure 12(a) and (b) are similar to Figure 11(a) and (b). Since Ay cT i > y c dev, the 
E/*i2f transition does not leave the violet area and the red area significantly increases in (c). Then the 
U* transition starts and the p distribution shifts to the right as shown in (d). The red area significantly 
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decreases from (c) to (d), and tension decreases and thus the "reversal of rate of tension recovery" 
cited at the beginning of this section is expected. The filament sliding continues in Phase 4, and the p 
distribution asymptotically approaches the isometric tetanus state (e). 

Figure 12. Changes of the p distribution during the isometric tension transient in the case 
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Figure 13 illustrates changes in Phases T1-T4 in another way by the 777o 0 bs vs. ALh s relation. The 
solid black arrows "Phase Tn" indicate an example of change of relative tension in Phase Tn for the case 
of Ay < ycdev and the dashed black arrows for the case of Ay > y c d ev . The edge Ay cT i in Figures 1 1(a) or 
12(a) is resulted from the elastic changes of the crossbridge and myosin filament, and is given as a 
function of ALh s by combining Equations 18 and 22: 

Ay cT i = -AL hs /2.2 (31) 

The tension variation Tri/Tbobs in Phase Tl can be calculated as a function of ALh s by using 
Equations 20 and 31. Calculation results are shown by the green curve in Figure 13. The edge y C T2 in 
Figures 1 1(b) or 12(b) are the same quantity as Ay cV i in Figure 10(c), since they result from the elastic 
changes of the crossbridge, myosin filament and actin filament. Then, from Equations 19 and 23, 
we have 



y cT 2 =-AL hs /4.6 



(32) 



Ay cT 2 = Ay cV i means TWToobs = 7vi/7b 0 b s - By using Equation 32 and replacing Tyi by Tj 2 and Ay cV i by 
Ay C T2 in Equation 21, 7T2/7o 0 b s can be calculated as a function of AL^ S . The calculation result is given 
by the brown curve in Figure 13. In the steady filament sliding, there is the definite relation between 
the tension T and the parameter y c as shown in Figure 8, where the red curve shows the relation for the 
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Hill-type filament sliding and the blue curve shows the relation when the filament sliding deviates 
from the Hill-type. The parameter y c in Figure 8 is related with Ay c in Figures 1 1 and 12 by the relation 
y c = JcOobs + Ay c since Ay c is the change of y c from the value at the isometric tension (y c o 0 bs)- In the 
steady filament sliding, all elastic elongations of the crossbridge, the myosin filament and the actin 
filament contribute to the muscle elongation and thus Ay c (=y c - y c o 0 b s ) is equal to -ALh s /4.6 as in 
Equation 32. The relative tension T/T 0 in Figure 8 can be converted to T/T 0obs by using the ratio 
ToobJTo = 0.88 (Equation 2). Based upon these considerations, the red and blue curves in Figure 8 are 
reproduced with the same colors in Figure 13. 



Figure 13. Phases in the isometric tension transient on the 777o 0 bs vs. ALh S relation. The 
thick black solid and dashed arrows, respectively, indicate Phases in the cases of 
AycTi < ycdev and Ay cT i > y C dev Actually the thick solid and dashed arrows are on the same 
vertical lines, respectively, except for the case of Phase Tl. 
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In Figure 13, the Phases are indicated by the solid black arrows "Phase Tn" in the case of Ay < y c d ev 
(the case of Figure 11), and by the dashed black arrows in the case of Ay > y c d ev (the case of Figure 12). 
The rapid tension change in Phase Tl occurs along the green curve as indicated by the solid or dashed 
arrows "Phase Tl". Phase T2 is a result of the elastic change of the actin filament and thus the arrows 
"Phase T2" start from the green curve and end on the brown curve. The main part of Phase T3a 
corresponds to the rapid tension increase from Figure 11(b) to (c) or from Figure 12(b) to (c). It is 
difficult to determine where Phase T3a turns into Phase T3b in Figure 13. Trial calculations, however, 
show that fairly good agreement can be obtained by assuming that Phase T3a ends at the blue curve or 
the red curve. Thus, in the case of Ay < y C d ev , the solid arrow "Phase T3a" is depicted with its tip on the 
blue curve. As discussed above, it is plausible that the red area in Figure 1 1(d) is nearly equal to that in 
(c). Taking this point into account, the solid arrow "Phase T3b" is depicted almost parallel to the 
abscissa, in accordance with the observation of the "extreme reduction of rate of tension recovery" 
noted at the beginning of this section. In the case of Ay > y c dev, as discussed above, it is plausible that 
the red area significantly decreases from Figure 12(c) to (d). Accordingly, the dashed arrow "Phase 
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T3b" is depicted downward, in accordance with the observation of the "reversal of rate of tension 
recovery " noted at the beginning of this section. After Phase T3b, the tension T approaches 7b 0 bs by 
the filament sliding as shown by the arrows "Phase T4". The above argument on the solid and dashed 
arrows "Phase T3b" suggests that the reversal of rate of tension recovery occurs when |ALh s | is large. 
Supporting this conclusion, Figure 1 of the paper by Julian and Sollins [19] shows that the reversal of 
the rate occurs when |AL hs | is large. 

Ford et al. [16] reported experimental data on variation of L/Loobs during 0-9 ms, as cited by green 
circles in Figure 14. The range of ALh s was +1.5 — 6.0 nm in their experiment. Since our model is not 
applicable for positive ALh s , the case of +1.5 nm is omitted in Figure 14. Now we shall try to reproduce 
these results by calculation. 

Figure 14. L/Toobs as functions of time t for various length change steps ALhs in the 
isometric tension transient. The origin of t is set as the moment that Phase Tl finishes. Green 
circles: Experimental data cited from Figure 23 of Ford et al. [16]. Red lines: Values of 
Li23a /Lbobs calculated by using Equations 34, 35, 36 and the parameter values in Equation 38. 
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In Figure 14, values of L/Loobs quickly increase within about 1 ms and then gradually approach to 
the value at 9 ms. The mean time duration needed for one U* transition is longer than 10 ms for 
TITq > 0.5 in Figure 5. Accordingly, effect of U* transition is not considered. The experimental length 
changes | ALh s I are smaller than 7.0 nm as seen in Figure 14. Hence, the following discussion is 
done referring to the solid arrows which are in this | ALh s I range. The tension changes in Phases Tl, 
T2 and T3a are denoted as T\, T 2 and L 3a , respectively. Also symbols T\ 2 and Li 2 3a are defined 
as follows: 
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T n = T y + T 2 (33) 

7123a = 7-12 + T 3a (34) 

According to the scheme in Figure 13, the change of T\ is given by the green curve, and T 2 changes 
from the green curve to the brown curve. In Figure 14, the origin of time t is set equal to 0 when the 
change of Phase Tl finishes. The parameters Cm2 and C f in 2 are defined as the initial and final values 
of T 2 /T 0obs , which are given by the green and brown curves, respectively in Figure 13. If T 2 is 
approximately expressed by a single decay constant Tn, TWToobs is given by, 

TVToobs = Cf m 2 - (C f ln 2 ~ C; n i2)exp(-?/ 7"T2) (35) 

The constant ri2 is related with the elastic change of the actin filament and is independent of ALh s . The 
solid arrow "Phase T3b" is drawn between the brown curve and the blue curve in Figure 13, so that the 
magnitude of tension of 777o 0 bs in T3b is given by the difference between these curves. It is uncertain 
when the U*n transition and thus Phase T3a start. To make calculation simple, it is assumed that it 
starts at the same moment as Phase T2, i.e., at t = 0. Also the change of T 3a is approximately expressed 
with a single decay constant 7r 3a : 

TWToobs = (C ful 3a ~ Cini3a){ 1 ~ exp(-?/r T3 a)} (36) 

where (Cini3a ~ CfMa) represents the magnitude of the tension variation. C; n i3 a and Cf in 3a are given by 
the brown and blue curves, respectively in Figure 13. 

Now the problem is how the parameter r T 3 a changes with AL hs . In the discussion on the 
force-velocity relation in [2], the average time t c for a myosin head at y c to cross over the potential 
barrier U*(y c ) is expressed as t c = (1/A) exp (U*(y c )/kT) (Equations 4-2-1 in [2]), and U*(y c ) is 
expressed by U*(y c ) = U*q - by c (Equations 4-2-2 in [2]). These formulae give t c = B exp (~cy c ) where 
B and c = b/kT are constants. In Figure 11(b) the boundary between the violet and red areas is 
indicated as —L + v C T2 + xq. The relative relation between this boundary and Al is similar to that 
between y c and A2. Then, in analogy to the above relation, with constants 5' and c\ the relation 
T3a= B' exp (—eye) is expected as an approximate expression, where y c ' = —L + y C T2+ xo (Figure 1 1(b)). 
Then, as L and x 0 are constants, r T 3 a in Equation 36 is approximately given by r T 3 a = 5exp(-c'v C T2), 
where c' is constant. Since v C T2 = ~ALh s /4.6 (Equation 32), this relation can be rewritten as, 

r T 3a = A 3a exp(c3 a AL hs ) (37) 

where A 3a and C3 a are constants. This relation implies that ri3a becomes small and the tension variation 
becomes fast when | ALh s I increases. 

Trial calculations were done for Ti23 a = Tn +r 3a (Equation 34) to explain the experimental 
observations by changing parameters tt 2 , A 3a and c 3a in Equations 35, 36 and 37. Fairly good 
agreement with the experimental data is obtained as shown in Figure 14 by using the following 
parameter values: 

7x2 = 0.7 ms (38a) 



A 3a = 3 ms 



(38b) 
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c 3a = 0.5 (l/nm) (38c) 

In Section 4, it is assumed that the elastic change of the actin filament occurs 0.4-1 ms after the length 
changes of the crossbridge and myosin filament. The time constant x T 2 = 0.7 ms is in consistency with 
this assumption. 

6. Isotonic Velocity Transient 

Isotonic velocity transients were studied by Civan and Podolsky [17], Huxley et al. [15] and Sugi 
and Tsuchiya [22,23]. A muscle was stimulated and initially held at a constant length. It was then 
released suddenly and allowed to shorten under a constant load. In this section, discussion will refer to 
the experimental data presented in Figure 3 of the article by Civan and Podolsky [17]. Responses of 
muscle to the sudden load change are classified into four Phases by Huxley [15]. Analogously, we use 
terminology "Phase Vn" as mentioned in Section 4. While Phase number of Huxley is related with 
time sequence of the length changes, Phase Vn is more related with molecular processes rather than 
the time sequence. Phase VI is the length change which simultaneously occurs with load change. 
Phase V2 is rapid early shortening. Phase V3 is extreme reduction or reversal of shortening speed. 
Phase V4 is responsible to the fact that the filament sliding velocity temporarily becomes larger than 
the steady value. Phases V2, V3 and V4 overlap each other in their time courses. 

Figure 15 shows changes of the p distribution in these four Phases. Figure 15 (a) is the p distribution 
just after Phase VI. Since y c = 0 at the isometric tension (cf. Figure 10 (a)), y c in Figure 15 (a) is equal to 
the variation from the isometric tension, Ay c vi (cf Figure 10 (c)). We shall discuss the case of 
Aycvi > Vcdev, since most experiments on the velocity transient were done for Ay c > y c dev 

In Figure 15, the vertical red arrows show time courses of Phases V2, V3 and V4. They overlap 
each other indicating the overlap of their time courses. The width of the violet area x 0 in Figure 15(a) 
is the same as in the isometric tension since the U*u transition does not occur yet. Then the U*u 
transition starts, i.e., Phase V2 starts. As mentioned above, there is a time delay in occurrence of the 
U* transition and thus of the filament sliding. Figure 15 (bl) shows the state in midway of Phase V2, 
where the U* transition starts. In (bl), a portion of the violet area has turned into red due to the U*uf 
transition and internal stress increases, which causes shrinkage of the p distribution, and thus of the 
muscle. At (b2), the U*m transition is over and all the violet area is turned into the red area as 
Ay c vi > Vcdev The process from (a) to (b2) is related with the U*u transition and is called Phase V2, 
where the red area increases and the internal stress increases. The width of p distribution changes from 
L in (a) to L-AL in (b2). Accordingly, a fast shortening of the muscle is expected. 

During Phase V2 the U* transition starts at (bl). There is a mutual interaction between the U* 
transition and the internal stress. The U* transition causes the filament sliding, releases the internal 
stress and tends to make y c larger. The internal stress tends to expand the p distribution, increases y c 
and accelerates the U* transition. The former process is called Phase V3. The latter process is called 
Phase V4. The mutual interaction almost disappears leaving a large value of y c at (c). Then the 
frequency of the U ^transition and y c change toward their steady values from (c) to (d), where the p 
distribution at the steady filament sliding is shown. Phases V3 and V4 occur in parallel starting with 
the U* transition at (bl) and end near (d). 
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Figure 15. Changes of the p distribution in the isotonic velocity transient when Ay c vi > y c dev- 
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Figure 16 illustrates changes in Phase V1-V4 in another way. This figure is drawn considering the 
case of small T/T 0o b s . (An example of the experimental data for small 777o 0 b s can be seen in 
Figure 18 (c) given later where AT/T 0o b s = 0.87, i.e., T/T 0o b s = 0.13.) The red, blue and brown curves 
are the same as the curves of the same colors in Figure 13. The black arrows indicate an example of 
changes of ALh S in the four Phases. As discussed above, there is the overlap between the Phases, and 
the arrow "Phase Vn" indicates the Phase which mainly contributes to the process. 

The length change in Phase VI occurs along the brown curve, since the elastic changes of the 
crossbridge, the myosin filament and the actin filament contribute to this change. The arrow "Phase 
V2" corresponds to the rapid decrease of muscle length from Figure 15 (a) to (b2). The length of the 
arrow "Phase V2" is tentatively depicted as 6 nm (the same order of magnitude of the rapid drop of 
5 nm in Figure 18 (c) given later). As discussed above, the time courses of Phases 3 and 4 overlap each 
other. The direction of arrow "Phase V3" is reversed, representing the muscle elongation from 
Figure 15 (b2) to (c). Through "Phase 4", the system approaches to the red curve corresponding to the 
change from (c) to (d) in Figure 15. The green arrow "Lhsv" symbolically represents the steady filament 
sliding corresponding to the state in Figure 15 (d). 

Now let us numerically reproduce the experimental data presented in Figure 3 of the paper by Civan 
and Podolsky [17]. The length change per half sarcomere at the steady filament sliding is denoted as Lh sv 
which corresponds to the dashed line in Figure 3 of [17]. Discussion will be done referring to L hsv 
expressed by 



Lhs 



-vt 



(39) 



Here v is the sliding velocity determined by the experimental data in [17]. 
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Figure 16. Phases in the isometric tension transient on the 777o 0 bs vs. ALh S relation. The 
black arrows indicate an example of changes of the relevant Phases in the case of 
Avcvi > Vcdev- Actually the arrows are on the same horizontal lines, except for the case of 
Phase VI. 
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The length change caused by the U*u transition in Phase V2 is denoted as Lh S 2- The speed of Lh S 2 
depends upon the frequency of U*n transition and will be large at the beginning and gradually decay. 
Its time course is approximately expressed with decay constant zv2 by 

L hs2 = -B V2 ( 1 - exp(-f/ r v2 )) (40) 

where By 2 is a constant. 

As discussed above, there is the mutual interaction between the U* transition and the internal stress, 
(a) The U* transition causes the filament sliding which releases the internal stress, (b) The internal 
stress tends to expand the p distribution, pushes y c forward and accelerates the U* transition. The 
interaction (a) is discussed first. The interaction (b) will be discussed later in relation with Phase 4. 

The interaction (a) releases the internal stress and thus elongates the muscle. As shown in Figure 15, 
there is an overlap between Phase V2 and V3, and thus there is an overlap between length change L^ s2 
and the length change due to the interaction (a). The combined length change is denoted as Lh S 23, 
which is expressed by multiplying Lh S 2 by exp(-?/ ZV3): 

£hs23 = -Bv2 { 1 - exp(-f/ r V 2)} exp(-?/ r V 3) (4 1 ) 

An example of L^ s2 3 is shown by the dotted black curve in Figure 17 to demonstrate its characteristics, 
together with Lh sv (the blue line). 

The decay times vy 2 and 7V3 should be functions of T/Toobs or v. The magnitude of Ty 2 is 
proportional to the time duration of occurrence of the U*n transition. Analogously to the manner to 
derive Equation 37, the relationx T3a = 5exp (— cy c ') is expected as an approximate expression, where y c ' 
is the y at the boundary between the violet area and red area in Figure 15(a), where y c ' = — L + Ay c vi+xo. 
The tension T becomes smaller as A y c vi becomes larger. As a simple approximation, AT = Tbobs ~ T is 
set proportional to A y cV i. Then the decay time Ty 2 is given approximately by 

7v2 = a V 2exp(c V 27"/7oobs) (42) 
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where av2 and Cv2 are constants. The time constant zv3 in Phase V3 is related with the U* transition, 
and its magnitude seems to be an order of r step = L/v (Equation 5). Hence r V 3 is set as 

r V3 = b w3 (L/v) (43) 

where Z?v3 is a constant. 

To see the characteristics of combination of Lh sv and Lh S 23, Lh S 23v is defined by 

£hs23v = £hsv + ^hs23 (44) 

An example of Lh S 23v is illustrated by the dashed violet curve in Figure 17. 

Figure 17. Characteristics of Lh s4 , Lh sv , Lhs23, £h S 23v and ALh s ' are illustrated in the case of 
Ar/r 0o bs = (T'oobs - T)/T 0o b s = 0.87, by using the equations and parameter values given in text. 
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Now let us consider about the part (b) of the interaction that the internal stress tends to expand the p 
distribution, pushes y c forward and accelerates the U* transition. This effect causes muscle elongation. 
In Figure 15, this effect is illustrated by depicting v cVc in (c) larger than y cV d in (d). In Phase 4, v cVc 
changes into v c vd and the filament sliding approaches to the steady value. As shown by the overlap of 
the arrows of Phases V3 and V4 in Figure 15, these effects overlap with each other. The length change 
in Phase V4 is denoted as L hs4 . In analogy to the expression of L hs 23 (Equation (41)), L hs 4 is 
approximately expressed by 

L hs 4 = 5v4 { 1 - exp(-?/ zv4) } 2 exp(-f/ r V s) (45) 

The characteristics of Lh S 4 expressed by this equation are illustrated by the dashed green curve in 
Figure 17. The internal stress is zero at t = 0 and increases by the term {1 - exp(-?/zv4)} 2 . The term 
Qxp(~t/ TVs) corresponds to the decrease of the muscle length due to the filament sliding. Concerning 
the magnitude of 5v4, it is plausible that the internal stress pushes y c forward more effectively when 
the internal stress rapidly increases, i.e., when the frequency of the U*n transition rapidly increases. 
The frequency of the U*n transition is given by the reciprocal of the duration rv2 given by 
Equation 42. Thus B y4 is expressed by 

5v4 = GW ZV2 (46) 



Int. J. Mol. Set 2011, 12 



1721 



where G V 4 is a constant. The time constants ZV4 and ZV5 will be mainly related with the U* transitions, 
and in analogy to Equation 43 they are expressed as 

r V 4 = b V 4(£/v) (47) 

r V5 = bv5(L/v) (48) 

where b* and bs are constant. 

Since the sum of L hsv , L hs 23 and L hs 4 is the total length change, it is denoted as AL hs ' 

ALhs' = Lh s 23 + Lhs4 + ^hsv (49) 

Note that ALh s ' is the length change from the moment when Phase 1 finishes, while the abscissa ALh S in 
Figure 16 includes the length change in Phase 1. 

Time course of ALh s ' was calculated with various trial values of the parameters looking for good 
agreement with the experimental data in the cases of (7b 0 bs - 7)/7o 0 b s = 0.22, 0.44 and 0.87. Calculation 
results with the following parameter values are shown in Figure 18. 

Figure 18. Length change ALh s ' vs. time t for three values of Ar/r 0o bs = (Toobs ~ T)/T 0o \, s . 
(a) AT/T 0ohs = 0.22; (b) AT/T^bs = 0.44; (c) AT/T 0ohs = 0.87. Green data points are obtained 
from Figure 3 of the paper by Civan and Podolsky [17]. Blue straight lines are Lh sv = ~vt 
(Equation 39). Red curves are AL^ calculated by using Equation 49 and the parameter 
values in Equation 50. 
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b Y4 = 0.8 (50f) 
bvs = 0.3 (50g) 

Characteristic features of the series of experimental data are fairly well reproduced by the 
calculation. The drastic change of the curve shape for different AT/T 0o b s is mainly due to the T 
dependence of the frequency of the occurrence of the U*n transition, i.e., the relation 
Tw2 = a V 2exp(c V 2777oobs) (Equation 42). 



7. Deviation from Hill-Type Force-Velocity Relation 

Let us consider the deviation from Hill-type force-velocity relation in connection of the 
£/*i2 transition. The deviation is determined by change of crossbridge population in well A2 and the 
U* potential barrier height. The maximum ratio of crossbridge population in Al per that in A2 is given 
by xo/(L-xo) in Figure 7 (b3). This ratio is not very large and the effect of population change is 
neglected in the following approximate formulation. 

The relative tensions T^/Ib and T*/To are given as functions of v c H and y c * in Figure 8. The same 
relationship can be expressed by giving y c H and y c * as functions of T/T 0 . From Equation 8, we have 

y c H = L[K t -{/a 2 - (Kf- Kb){K t - Ipof'KToL)} m }]/(/Cf- K h ) (51) 

From Equation 12, we have 

y * = [{b-(T*/T 0 )}/a] m (52) 

These relations and y c u - y c * are shown as functions of T/T 0 = T^/Tq = T*/T 0 in the range between 
Tdev/To and 7oobs/7o in Figure 19 (a). 

The filament sliding velocity of Hill-type is denoted as v H and the velocity in the presence of the 
C/*i2 transition as v*. The period for the crossbridge to move over L is given as r ste p = L/v u in Figure 5, 
which should have a similar nature as the time duration of U* transition, t c = (l/A)exp(t/*/&7) 
(Equation 4-2-1 in [2]), where A is a constant and U* is the potential barrier for the U* transition. Thus 
v H = L/r ste p is approximately proportional to l/exp(£/*//c7). Here U* is the potential barrier height when 
there is no U*n transition. The difference y c H - y c * shown in Figure 19 (a) can be regarded as a 
measure of the magnitude of structural change in MA 3 . The potential barrier U* should be affected by 
this structural change, and the change of the U* height is approximately set as bv(y c H ~ y c *) where by 
is a constant proportional to 1/kT. Then Qxp(U*/kT) changes into exp((U*/kT) + (by(y c u - y c *)) in 
Tdev <T< Toobs- Then, if we put v H = C/exp(U*/kT) and v* = C/exp((U*/kT) + K(y c n - y 0 *)) with a 
constant C, we have 

v* = v H /exp(Mvc H -v c *)) (53) 

Since v c H = y c * = v c d ev (cf. Figure 8), Equation 53 gives v* = v H at v c dev, as required. The blue curve in 
Figure 19(b) shows results of calculation by Equation 53 with the parameter value 

fc v = 7.9(l/nm) (54) 



Int. J. Mol. Sci. 2011, 12 



1723 



Agreement with the experimental data is fairly good. Accordingly, the deviation from Hill-type 
force-velocity relation is mainly due to the change of the potential barrier height U* caused by 
occurrence of U*u transition. 

Figure 19. Explanation of the deviation of the force-velocity relation from Hill-type, 
(a) v c H , v c *, v c H - y c * as functions of 77Tb, calculated by Equations 51 and 52; (b) Brown 
circles: Experimental data of the sliding velocity per half sarcomere cited from Fig, 3A of the 
paper by Edman [12]. Red curve: v H , the Hill-type velocity calculated in [2]. Blue curve: v*, 
the velocity in the presence of the U*n transition calculated by Equations 53 and 54. 
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8. Summary and Discussion 

In our previous paper [2] (the first part of this Remarks series), difficulty of the power stroke model 
is pointed out and a new model is proposed to avoid the difficulty. In the model, it is proposed that 
about 41% of the myosin heads are bound to actin filament and each bound head forms complex MA3 
with three actin molecules. The complex MA 3 translates along the actin filament changing its partner 
actin molecules in cooperation with U* transition. This model well explains the properties in the 
steady filament sliding such as the tension-dependence of the muscle stiffness, the Hill-type force 
velocity relation and the tension-dependence of energy liberation rate. In the present paper, the 
isometric tension transient and isometric velocity transient are studied based upon the model. 
Statistical ensemble of crossbridges is considered and the binding probability density p is introduced. 
On rectangular p approximation, the edge of the rectangle, y c determines dynamic properties of muscle 
in cooperation with U* and U*n transitions. The internal structure of the MA3 complex becomes 
temporally unstable by the sudden length change or by the sudden load change. The complex muscle 
behaviors observed in these transients are related with the process that the disturbed internal structure 
returns to its stationary state. 
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Results reported in the present paper are summarized as follows. 

(1) The tension variations in the first Phases in the isometric tension transient (Ford et al. [16]) 
and the isotonic velocity transient (Civan and Podolsky [17]) are well explained as shown 
in Figure 9. 

(2) Ratios of extensibilities of crossbridge, myosin filament and actin filament are estimated as 
0.22, 0.26 and 0.52 (Equation 30), in reasonable agreement with the approximate values 
(0.31, 0.27, 0.42) determined by X-ray diffraction by Wakabayashi et al. [21]. 

(3) The experimental data on the isometric tension transient reported by Ford et al. [16] are fairly 
well explained as shown in Figure 14. 

(4) The characteristic features of muscle in the isotonic velocity transient observed by Civan and 
Podolsky [8] are fairly well explained as shown in Figure 18. 

(5) The deviation from the Hill-type force-velocity relation observed by Edman [12] is reproduced 
as shown in Figure 19(b). 

It should be noted that the above-mentioned agreements between experimental data and calculation 
results are obtained by using the muscle stiffnesses k? and k\, determined in [2], whose numerical 
values are given in Appendix. 

The obtained results suggest that the ideas of the ensemble of crossbridges and the rectangular 
p approximation are useful tools in theoretical studies of muscle contraction. 

In Figure 8, tangent of the T*ITq vs. y c * relation (the blue curve) is 0 at y c * = 0. Accordingly, small 
fluctuation of tension around r 0o bs can produce significant variation of the muscle length. In this 
connection, the spontaneous oscillatory contraction (SPOC) of muscle (cf. the review by Ishiwata and 
Yasuda [24]) seems interesting. As mentioned in [24], it is a possibility that SPOC has some relation 
with the activities of cardiac muscles (cf. [25]). 
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Appendix. Symbols and Their Values 

kT: 3.77 X10" 21 atO°C. 

Kf. stiffness of an elongated crossbridge. 2.80 pN/nm (Equation 4-1-13 in [2]). 

Kb: stiffness of a shrinked crossbridge. 0.26 pN/nm (Equation 4-1-14 in [2]). 

L: period of binding sites of the myosin head along the actin filament, 5.46 nm. 

Lh s : length of the half sarcomere. 

p: mean force produced by one crossbridge. 

po: maximum p in the extended Hill-type force-velocity relation. 

Poobs: P in isometric contraction. 5.7 pN (Equation 2-2-3 in [2]). 

r. ratio = (number of myosin heads attached to actin filaments)/(total number of myosin heads) in 
contracting muscle. 0.41 (Equation 3-1-1 in [2]). 
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tcb: extensibility ratio of crossbridge. 0.22 (Equation 30a). 

r M : extensibility ratio of myosin filament. 026 (Equation 30b). 

r A : extensibility ratio of actin filament. 0.52 (Equation 30c). 

s: sarcomere length of muscle with a full filament overlap. 2.10 um [3]. 

T: tension (P is used for T in [2]). 

T 0 : maximum T in the extended Hill-type force-velocity relation (P 0 in [2]). 

5 2 

7oob s : isometric tension of muscle with a full filament overlap. 4.1 x 10 N/m [3]. 7o 0 b s /7o = 0.88 
(cf. Figure 1). 

Td ev : T at which the velocity starts to deviate from Hill-type force-velocity relation. 

7dev/7o = 0.66 (cf. Figure 1). 

T* \ observed tension for T> 7d ev . 

Vmax'. velocity of the filament sliding under no load on muscle. 2.36 um/s at 1.8 °C [3]. 
y: srinkage shortening of the crossbridge (cf. Figure 4). 
v c : maximum y in the rectangular p distribution. 
v c (0): y c for v = v max . 4.2 nm (Equation 4-1-11 in [2]). 

Vco: v c for v = 0 in the Hill-type force-velocity relation. 0.73 nm (Equation 4-1-12 in [2]). 

Vcdev: Vc for T dev /T 0 = 0.66. 1.60 nm (Figure 8). 

y c oob s : y c for T 0ohs /T 0 = 0.88. 1.03 nm (Figure 8). 

Vc*: v c observed in the presence of the U*u transition (Figure 7). 

References 

1 . Mitsui, T. Induced potential model of muscular contraction mechanism and myosin molecular 
structure. Adv. Biophys. 1999, 36, 107-158. 

2. Mitsui, T.; Ohshima, H. Remarks on muscle contraction mechanism. Int. J. Mol. Sci. 2008, 9, 
872-904. 

3. Oshima, K.; Takezawa, Y.; Sugimoto,Y.; Kobayashi, T.; Irving, T.C.; Wakabayashi, K. Axial 
dispositions and conformations of myosin crossbridges along thick filaments in relaxed and 
contracting states of vertebrate striated muscles by X-ray fiber diffraction. J. Mol. Biol. 2007, 367, 
275-301. 

4. Knupp, C; Offer, G.; Ranatunga, K.W.; Squire, J.M. Proving muscle myosin motor action: X-ray 
(M3 and M6) interference measurements report motor domain not lever arm movement. J. Mol. 
Biol. 2009, 390, 168-181. 

5. Andreeva, A.L.; Andreev, OA.; Borejdo, J. Structure of the 265-kilodalton complex formed upon 
EDC cross-linking of subfragment 1 to F-actin. Biochemistry 1993, 32, 13956-13960. 

6. Podolsky, R.J.; Onge, S.S.; Yu, L. X-ray diffraction of actively shortening muscle. Proc. Natl. 
Acad. Sci. USA 1976, 73, 813-817. 

7. Huxley, H.E. Time resolved X-ray diffraction studies in muscle. In Cross-bridge Mechanism in 
Muscle Contraction; Sugi, H., Pollack, G.H., Eds.; University of Tokyo Press: Tokyo, Japan, 
1979;pp.391-405. 

8. Huxley, H.E.; Kress, M. Crossbridge behaviour during muscle contraction. J. Muse. Res. Cell 
Motility 1985, 6, 153-161. 



Int. J. Mol. Sci. 2011, 12 



1726 



9. Yagi, N.; Takemori, S.; Watanabe, M. An X-ray diffraction study of frog skeletal muscle during 
shortening near the maximum velocity. J. Mol. Biol. 1993, 231, 668-677. 

10. Hill, A.V. The heat of shortening and the dynamic constants of muscle. Proc. Roy. Soc. 1938, 
B126, 136-195. 

11. Hill, A.V. The effect of load on the heat of shortening muscles. Proc. Roy. Soc. B 1964, 159, 
297-318. 

12. Edman, K.A.P. Double-hyperbolic force-velocity relation in frog muscle fibres. J. Physiol. 1988, 
404,301-321. 

13. Eyring, H. The Eyring equations. In Physical Chemistry, 6th ed.; Atkins, P.W., Ed.; Sumanas, 
Inc.: Pasadena, CA, USA, 1998. 

14. Huxley, A.F.; Simmons, R.M. Proposed mechanism of force generation in striated muscle. Nature 
1971, 233, 533-538. 

15. Huxley, A.F. Muscular contraction. J. Physiol. 1974, 243, 1-43. 

16. Ford, L.E.; Huxley, A.F.; Simmons, R.M. Tension responses to sudden length change in 
stimulated frog muscle fibres near slack length J. Physiol. 1977, 269, 441-515. 

17. Civan, M.M.; Podolsky, R.J. Contraction kinetics of striated muscle fibres following quick 
changes in load. J. Physiol. 1966, 184, 511-534. 

18. Irving, M. Give in the filament. Nature 1955, 374, 14-15. 

19. Julian, F.J.; Sollins, M.R. Variation of muscle stiffness with force at increasing speeds of 
shortening. J. Gen. Pysiol. 1975, 66, 287-302. 

20. Huxley, H.E.; Stewart, A.; Sosa, H.; Irving, T. X-ray diffraction measurements of the extensibility 
of actin and myosin ffilaments in contracting muscle. Biophys. J. 1994, 67, 241 1-2421. 

21. Wakabayashi, K.; Sugimoto, Y.; Tanaka, H.; Ueno, Y.; Takazawa, Y.; Amemiya, Y. X-ray 
diffraction evidence for the extensibility of actin and myosin filaments during muscle contraction. 
Biophys. J. 1994, 67, 2422-2435. 

22. Sugi, H.; Tsuchiya, T. Isotonic velocity transients in frog muscle fibres following quick changes 
in load. J. Physiol. 1981, 319, 219-238. 

23. Sugi, H.; Tsuchiya, T. Enhancement of mechanical performance in frog muscle fibres after quick 
increases in load. J. Physiol. 1981, 319, 239-252. 

24. Ishiwata, S.; Yasuda, K. Mechano-chemical coupling in spontaneous oscillatory contraction of 
muscle. Phase Trans. 1993, 45, 105-136. 

25. Bers, D. Exitation contraction coupling and cardiac contractile force. In Development in 
Cardiovascular Medicine, 2nd ed.; Springer: Berlin, Germany, 2001; Volume 237. 

© 2011 by the authors; licensee MDPI, Basel, Switzerland. This article is an open access article 
distributed under the terms and conditions of the Creative Commons Attribution license 
(http://creativecommons.Org/licenses/by/3 .0/) . 



